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Abstract 


A  vectorial  Fourier-based  technique  for  electromagnetic  (EM)  wave  reconstruction  with 
application  to  polarimetric  airborne  and  spacebome  radar  data  exploitation  is  presented.  The 
approach  provides  a  comprehensive  picture  of  scattering  through  an  electromagnetic-based 
formulation  of  equivalent  sources  and  introduction  of  vector  scattering  functions.  The  method 
is  different  from  conventional  modeling  techniques  for  SAR  applications  as  result  of  the  full 
electromagnetic  treatment  of  field  interaction  with  the  scatterer,  the  possibility  of  introducing 
new  and  controllable  feature  classes  for  target  classification,  and  the  accurate  decomposition 
of  source  impulse  response  function  that  avoids  potential  errors  (e.g.,  loss  of  coherent 
information)  associated  with  the  spherical  phase  approximations.  In  addition  to  including  the 
potential  for  high  quality  image  formation,  the  capability  of  extracting  target  or  scatterer 
information,  such  as  coherent  radar  cross  section  (RCS)  is  explored.  It  is  highlighted  how 
target  image  and  information  processing  can  be  optimized  and  tailored  for  a  desired  scenario 
(e.g.,  the  required  information  type,  the  nature  of  the  targets)  by  applying  the  developed 
model  and  devising  appropriate  inversion  techniques  to  aid  data  processing  in  support  of 
exploitation. 


Resume 


Le  document  presente  une  technique  vectorielle  a  transformee  de  Fourier  permettant  de 
reconstituer  des  ondes  electromagnetiques  (EM)  par  l'application  a  I'exploitation  des  donnees 
de  radars  polarimetriques  aeroportes  et  spatiaux.  Cette  technique  procure  une  image  globale 
de  la  diffusion  par  la  formulation  de  sources  equivalentes,  selon  des  principes 
electromagnetiques,  et  le  recours  a  des  fonctions  vectorielles  de  diffusion.  Elle  differe  des 
techniques  classiques  de  moderation  des  applications  SAR  en  raison  du  traitement 
electromagnetique  complet  de  l'interaction  des  champs  avec  le  diffuseur,  de  la  possibilite 
d'etablir  des  ensembles  de  caracteristiques  nouveaux  et  controlables  pour  la  classification  des 
cibles  et  de  la  decomposition  precise  de  la  fonction  de  reponse  impulsive  des  sources,  qui 
elimine  les  erreurs  possibles  (p.  ex.  la  perte  d'information  coherente)  liees  aux  approximations 
de  phase  spherique.  Outre  la  formation  d'images  de  grande  qualite,  la  possibilite  d'extraire  de 
I'information  sur  les  cibles  ou  les  diffuseurs,  comme  la  surface  equivalente  radar  (SER) 
coherente,  est  egalement  examinee.  Le  document  fait  ressortir  comment  le  traitement  des 
images  des  cibles  et  de  I'information  sur  les  cibles  peut  etre  optimise  et  adapte  en  fonction 
d'un  scenario  desire  (p.  ex.,  le  type  de  I'information  requise  et  la  nature  des  cibles)  par 
l'application  du  modele  elabore  et  la  conception  de  techniques  d'inversion  appropriees 
facilitant  Ie  traitement  des  donnees  a  l'appui  de  I'exploitation. 
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Executive  summary 


Polarimetric  radar  or  SAR  (Synthetic  Aperture  Radar)  detection  is  a  powerful  remote  sensing 
technique  for  potential  military  applications  concerning  target  recognition,  identification  and 
classification.  The  polarimetric  approach  is  potent  as  it  effectively  explores  radar  target 
features  (e.g.,  structure,  change)  by  making  use  of  the  full  vector  properties  of  the 
electromagnetic  (EM)  field  while  interacting  with  the  target  scatterers.  Accordingly, 
exploitation  of  polarimetric  radar  data,  using  appropriate  algorithms,  results  in  an  enhanced 
target  detection  and  classification  capability  compared  to  single  channel,  single  polarization 
data  exploitation.  Such  improved  functionality  is  desired  for  consistent  target  detection, 
reduction  of  misclassification  probability,  reduction  of  false  alarm  rate  (FAR),  and  reliable 
target  recognition. 

It  is  evident  from  the  characteristic  description  of  radar  polarimetry  that  the  analysis  of  target 
scattering  ought  to  be  based  on  electromagnetic  theory.  A  full  electromagnetic  approach  to 
modeling  the  scattered  vector  fields  and  characterize  the  targets  introduces  new  challenges 
and  requirements  for  more  rigorous  algorithm  development.  On  the  other  hand,  the  theoretical 
framework  cannot  be  very  complicated  from  a  practical  application  point  of  view.  The 
optimal  trade-off  between  these  two  aspects,  i.e.,  accuracy  and  feasibility,  is  an  important 
issue  that  has  not  been  systematically  addressed  in  the  literature  and  is  the  underlying 
motivation  for  this  work. 

Considering  the  diversity  of  features  characterizing  a  radar  or  SAR  scenario  for  remote 
sensing  (e.g.,  application  types,  required  analysis  types,  target  and  clutter  types),  development 
of  a  unique  optimized  approach  to  model  the  scattering  for  polarimetric  applications  is  a  huge 
if  not  insurmountable  task.  The  intention  is  not  to  accomplish  such  here  as  this  is  not  realistic. 
Instead,  the  objective  is  to  construct  a  theoretical  framework  capable  of  accurately  modeling 
the  process  of  vector  EM  field  interaction  with  the  target  scatterer.  This  general  formulation 
will  be  tailored  based  on  the  requirements  for  the  application  classes  of  interest  and  will  be 
applied  for  simulations  and  experimental  purposes  in  later  works.  The  contributions  of  this 
work  can  be  viewed  in  two  folds: 

The  first  contribution  is  a  thorough  survey  of  the  literature  in  polarimetric  remote  sensing  for 
modeling  of  EM  scattering  and  characterizing  the  target  scatterer.  The  outcome  signifies 
nature,  objectives,  required  deliverables,  major  classes  and  capabilities/limitations  of  the 
existing  analytic  techniques  for  radar  target  analysis.  In  parallel,  a  study  is  performed  to 
identify  robust  and  efficient  methods  to  model  scattering  and  diffraction  (as  an  EM 
phenomenon)  that  belong  to  different  areas  within  the  electromagnetic  realm  (e.g.,  microwave 
structures)  but  have  similar  theoretical  foundations.  The  overall  objective  from  the  above 
investigations  is  to  derive  a  set  of  key  requirements  for  the  modeling  formulation,  and  at  the 
same  time,  propose  a  feasible  implementation  of  the  associated  algorithms. 

The  second  contribution  of  this  work  is  the  development  of  a  theoretical  approach  to  model 
the  target  scattering  (based  on  the  key  requirements  identified)  that  is  different  from 
conventional  methods  in  two  ways.  One  aspect  is  the  rigorous  analysis  (usually  lacking  in 
typical  analyses)  and  application  of  source  Green’s  function  through  a  four-dimensional  (i.e., 
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space-time)  Fourier  decomposition.  When  applying  a  coherent  approach  for  target 
identification,  this  accurate  spectral  treatment  of  Green’s  function  with  respect  to  both  phase 
and  amplitude  provides  the  means  for  more  accurate  modeling  of  EM  field  scattering  by  the 
target  scatterer.  The  other  aspect  is  the  representation  and/or  classification  of  target  scatterers 
by  vector  current  functions.  This  approach,  which  is  based  on  the  concept  of  equivalent 
sources  in  electromagnetic  theory,  is  close  to  the  roots  of  EM  scattering  governed  by 
Maxwell’s  equations.  Technically,  the  latter  is  on  the  same  basis  as  the  conventional 
scattering  matrix  analysis  due  to  formulation  of  scattering-induced  sources.  However,  the 
present  approach  is  novel  (to  the  best  of  author’s  knowledge  in  the  field  of  polarimetric  SAR 
scattering  modeling  for  target  classification)  and  appealing  from  the  following  points  of  view. 

To  classify  the  target  scatterer,  the  associated  vector  current  function  in  the  space  or  spectral 
domain  is  decomposed  into  a  set  of  basis  functions  that  can  be  viewed  as  features.  It  is 
possible  to  perform  this  decomposition  using  different  basis  sets  where  each  set  represents  a 
class  of  features  or  sub-features.  Considering  the  generic  nature  of  problem  formulation, 
efficient  techniques  may  be  adopted  from  other  areas  (e.g.,  RF,  microwaves)  for  such  target 
characteristic  function  decompositions.  As  can  be  seen,  there  exist  degrees  of  freedom  and  the 
possibility  of  defining  various  basis  sets  and  number  of  members  or  basis  functions  for  each 
set.  The  described  potential  of  controlling  the  nature  and  number  of  features  for  target 
decomposition  based  on  the  application  requirements,  alleviates  the  misclassification 
problems  associated  with  coarse  classification/segmentation  encountered  through  current 
coherent  target  decomposition  techniques.  It  is  also  conceivable  that  study  of  the  induced 
equivalent  current  source’s  behavior  on  the  scatterer,  yield  a  more  consistent  and  tangible 
understanding  of  the  scatterer’s  physics.  Finally,  the  source  functions  derived  through  the 
present  approach  may  be  readily  applied  to  calculate  the  associated  EM  field  quantities  (such 
as  coherent  radar  cross  section)  in  both  near  and  far  field  zones. 

The  present  methodology  provides  a  core  theoretical  framework  for  polarimetric  SAR 
modeling.  The  approach  is  general  and  should  be  tailored  by  investigating  the  SAR  scenario 
of  interest  and  designing  a  suitable  structure  for  required  feature  classes.  A  significant 
advantage  is  the  possibility  of  introducing  various  target  classification  or  basis  sets  to  the  core 
formulation.  Moreover,  is  the  opportunity  of  adopting  robust  and  well-proven  techniques 
(e.g.,  method  of  moments,  image  theory)  from  other  EM  areas  to  obtain  the  optimum  (i.e., 
best  target  representative)  basis  sets  and  associated  decomposition  mechanisms.  The  future 
work  will  be  directed  towards  model  or  sub-model  development  for  different  types  of 
applications.  Various  EM-based  numerical  techniques  should  be  studied  and  properly  tailored 
for  desired  applications.  Convenient  experimental  schemes  must  still  be  designed  to  validate 
and  tune  the  model. 


Sabry,R.  2004.  Electromagnetic  modeling  in  the  spectral  domain  for  polarimetric  radar 
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Sommaire 


La  detection  par  radar  polarimetrique  ou  SAR  (radar  a  synthese  d'ouverture)  constitue  une 
puissante  technique  de  teledetection  pour  les  applications  militaires  possibles  liees  a  la 
reconnaissance,  a  l'identification  et  a  la  classification  des  cibles.  La  technique  polarimetrique 
s'avere  particulierement  efficace  du  fait  qu'elle  examine  avec  soin  les  caracteristiques  des 
cibles  radar  (comme  les  structures  et  les  changements)  en  tirant  parti  de  toutes  les  proprietes 
vectorielles  des  champs  electromagnetiques  (EM)  durant  leur  interaction  avec  les  diffuseurs 
de  cibles.  Sexploitation  des  donnees  de  radar  polarimetrique  au  moyen  d'algorithmes 
appropries  permet  done  d'ameliorer  la  detection  et  la  classification  des  cibles 
comparativement  a  l'exploitation  des  donnees  avec  canal  unique  et  polarisation  unique.  Cette 
amelioration  est  souhaitable  pour  la  detection  uniforme  des  cibles,  la  reduction  de  la 
probability  d'erreurs  de  classification,  la  reduction  du  taux  de  fausses  alarmes  (TFA)  et  la 
reconnaissance  fiable  des  cibles. 

D'apres  la  description  caracteristique  de  la  polarimetrie  radar,  il  est  evident  que  Panalyse  de  la 
diffusion  des  cibles  devrait  se  fonder  sur  la  theorie  electromagnetique.  Une  technique 
totalement  electromagnetique  de  moderation  des  champs  vectoriels  diffuses  et  de 
caracterisation  des  cibles  suscite  de  nouveaux  defis  et  impose  de  nouvelles  exigences  en  vue 
de  I'etablissement  d'algorithmes  plus  rigoureux.  Par  ailleurs,  le  cadre  theorique  ne  peut  pas 
etre  tres  complique  du  point  de  vue  de  l'application  pratique.  Le  compromis  optimal  entre  ces 
deux  aspects,  soit  la  precision  et  la  faisabilite,  represente  un  facteur  important  qui  n'a  pas 
encore  ete  traite  de  fa9on  systematique  dans  la  documentation  publiee  et  qui  constitue  la 
motivation  sous-jacente  a  ces  travaux. 

Etant  donne  la  diversity  des  caracteristiques  des  radars  ou  des  scenarios  SAR  de  teledetection 
(p.  ex.,  types  d'applications,  types  d'analyses  requises  et  types  de  cibles  et  de  clutter), 
I'etablissement  d'une  methode  optimisee  unique  permettant  de  modeliser  la  diffusion  des 
applications  polarimetriques  constitue  une  tache  immense,  voire  un  obstacle  insurmontable.  II 
serait  irrealiste  de  tenter  I'adoption  d'une  telle  methode.  L'objectif  consiste  plutot  a  etablir  un 
cadre  theorique  permettant  de  modeliser  avec  exactitude  le  processus  d'interaction  des  champs 
EM  vectoriels  avec  le  diffuseur  de  cible.  Cette  formulation  generate  sera  adaptee  aux 
exigences  des  classes  d'applications  d'interet  et  sera  appliquee  dans  des  travaux  ulterieurs  a 
des  fins  experimentales  et  de  simulation.  Les  contributions  de  ces  travaux  sont  doubles. 

La  premiere  contribution  prend  la  forme  d'un  examen  approfondi  de  la  documentation  en 
teledetection  polarimetrique  pour  la  modelisation  de  la  diffusion  EM  et  la  caracterisation  du 
diffuseur  de  cible.  Le  resultat  fait  ressortir  la  nature,  les  objectifs,  les  produits  livrables  requis 
ainsi  que  les  principales  classes  et  capacites/limites  des  techniques  analytiques  existantes  pour 
l'analyse  des  cibles  radar.  Parallelement,  une  etude  vise  a  determiner  des  methodes  robustes  et 
efficaces  de  modelisation  de  la  diffusion  et  de  la  diffraction  (en  tant  que  phenomenes  EM) 
ayant  trait  a  differents  aspects  de  I'electromagnetisme  (p.  ex.  les  structures  hyperfrequences) 
mais  partageant  des  fondements  theoriques  similaires.  Ces  recherches  ont  pour  objectif 
general  d'etablir  un  ensemble  d'exigences  cles  pour  la  formulation  de  la  modelisation  et, 
simultanement,  a  proposer  un  moyen  de  mise  en  oeuvre  des  algorithmes  connexes. 
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V 


A  titre  de  seconde  contribution,  ces  travaux  permettent  d'elaborer  une  methode  theorique 
servant  a  modeliser  la  diffusion  des  cibles  (d'apres  les  exigences  cles  etablies),  qui  differe  des 
methodes  habituelles  a  deux  egards.  En  premier  lieu,  cette  difference  tient  a  l'execution  d'une 
analyse  rigoureuse  (d'une  rigueur  superieure  a  celle  des  analyses  typiques)  et  a  1'application  de 
la  fonction  de  Green  pour  la  source,  par  une  decomposition  de  Fourier  quadridimensionnelle 
(espace-temps).  Si  Ton  recourt  a  une  methode  coherente  d'identification  des  cibles,  ce 
traitement  spectral  precis  de  la  fonction  de  Green  par  rapport  a  la  phase  et  a  l'amplitude 
permet  une  modelisation  plus  precise  de  la  diffusion  du  champ  EM  par  le  diffuseur  de  cible. 
En  second  lieu,  cette  methode  se  differencie  par  la  representation  et/ou  la  classification  des 
diffuseurs  de  cibles  au  moyen  de  fonctions  vectorielles  de  courant.  Une  telle  fa?on  de 
proceder,  fondee  sur  le  principe  des  sources  equivalentes  dans  la  theorie  de 
i'electromagnetisme,  se  rapproche  des  elements  de  la  diffusion  electromagnetique  regis  par  les 
equations  de  Maxwell.  Sur  le  plan  technique,  elle  se  trouve  sur  le  meme  pied  que  l'analyse 
matricielle  classique  de  la  diffusion,  a  cause  de  la  formulation  de  sources  induites  par  la 
diffusion.  La  methode  actuelle  est  toutefois  novatrice  (autant  qu'en  sache  l'auteur  dans  le 
domaine  de  la  modelisation  de  la  diffusion  par  SAR  polarimetrique  pour  la  classification  des 
cibles)  et  interessante  pour  les  raisons  decrites  ci-dessous. 

Pour  la  classification  des  diffuseurs  de  cibles,  la  fonction  vectorielle  de  courant  connexe  dans 
le  domaine  spatial  ou  spectral  est  decomposee  en  un  ensemble  de  fonctions  de  base  pouvant 
etre  considerees  comme  des  caracteristiques.  11  est  possible  d'effectuer  cette  decomposition  a 
l'aide  de  differents  ensembles  de  base,  chaque  ensemble  representant  une  classe  de 
caracteristiques  ou  de  sous-caracteristiques.  Compte  tenu  de  la  nature  gen6rique  de  la 
formulation  du  probleme,  il  est  possible  d'adopter  des  techniques  efficaces  issues  d'autres 
secteurs  (RF,  micro-ondes,  etc.)  pour  effectuer  de  telles  decompositions  de  fonctions 
caracteristiques  des  cibles.  Comme  on  peut  le  constater,  il  existe  des  degres  de  liberte  et  la 
possibility  de  definir  divers  ensembles  de  base  ainsi  qu'un  certain  nombre  de  membres  ou 
fonctions  de  base  pour  chaque  ensemble.  La  possibility  decrite  de  determiner  la  nature  et  le 
nombre  de  caracteristiques  de  decomposition  des  cibles  en  fonction  des  exigences  de 
I'appiication  attenue  les  problemes  de  classification  liys  a  la  classification/segmentation 
approximative,  qui  sont  typiques  des  techniques  actuelles  de  decomposition  coherente  des 
cibles.  On  peut  aussi  concevoir  que  1'etude  de  l'incidence  d'une  source  de  courant  equivalent 
induit  sur  le  diffuseur  procure  une  comprehension  plus  uniforme  et  tangible  de  la  physique  du 
diffuseur.  Enfin,  les  fonctions  de  source  liees  a  la  methode  actuelle  peuvent  s’appliquer 
directement  au  calcul  des  valeurs  des  champs  EM  connexes  (comme  la  surface  equivalente  de 
radar  coherent),  tant  en  champ  proche  qu'en  champ  lointain. 

La  presente  methode  fournit  un  cadre  theorique  central  pour  la  modelisation  par  SAR 
polarimetrique.  Cette  methode  est  generale  et  devrait  etre  adaptee  par  I’examen  du  scenario 
SAR  d'intyret  ainsi  que  par  la  conception  d'une  structure  appropriee  applicable  aux  classes  de 
caracteristiques  requises.  Un  avantage  significatif  tient  a  la  possibility  de  recourir  a  diverses 
classifications  de  cibles  ou  a  differents  ensembles  de  base  pour  la  formulation  centrale.  Il  est 
egalement  possible  d'adopter  des  techniques  robustes  ayant  fait  leurs  preuves  (methode  des 
moments,  theorie  des  images,  etc.)  dans  d'autres  secteurs  de  l'EM  afin  d'obtenir  les  ensembles 
de  base  optimaux  (le  plus  representatifs  des  cibles)  et  les  mecanismes  de  decomposition 
correspondants.  Les  travaux  futurs  viseront  I'etablissement  d'un  modele  ou  d'un  sous-modele 
pour  les  differents  types  d'applications.  Diverses  techniques  numeriques  fondees  sur 
i'electromagnetisme  devraient  etre  etudiees  et  adaptees  correctement  aux  applications 
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desirees.  II  reste  encore  a  concevoir  des  structures  experimentales  appropriees  pour  valider  et 
raffiner  le  modele. 
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1.  Introduction 


Radar  polarimetry  is  a  fast  growing  and  powerful  remote  sensing  technique.  It  is  anticipated 
that  the  majority  of  new  radars  will  have  polarimetric  or  scattering  matrix  capability  [1].  The 
strength  of  the  polarimetric  approach  is  due  to  many  important  and  key  features  of  radar 
targets  being  vectorial  in  nature  [2].  Hence,  making  use  of  the  inherent  vector  property  of 
electromagnetic  field  (EM)  interactions  with  the  scatterer,  results  in  a  more  comprehensive 
picture  of  the  target.  There  exists  a  considerable  volume  of  fundamental  work  in  classical 
physics  that  explores  different  aspects  of  electromagnetic  wave  scattering  and  diffraction  by 
various  objects  that  support  the  polarimetric  approach.  There  is  also  significant  research 
literature  on  polarimetric  applications  and  their  advantages  for  active  microwave  remote 
sensing  [3]-[6]. 

Polarimetric  detection  in  radar  technology  is  used  to  create  high-resolution  images  [7]-[9]  and 
provide  better  target  detection  capabilities  in  comparison  with  non-polarimetric  or  single 
channel  operation  [10]-[13].  The  latter,  which  amounts  to  the  analysis  of  a  detected  signal  to 
determine  the  EM  characteristic  (i.e.,  the  EM  signature)  of  the  target  for  target  identification, 
is  an  important  radar  application.  To  properly  correlate  the  received  EM  signal(s)  with  a 
particular  object  or  scatterer  (or  combination  of  scatterers),  an  EM-based  formulation  of 
scattering  and  propagation  should  be  applied. 

An  important  and  early  insight  to  radar  polarimetry  for  target  detection  [14]  was  to  represent  a 
target  as  a  polarization  transformer  and  to  express  its  coherent  properties  with  a  2x2 
scattering  matrix.  This  concept  was  developed  by  many  researchers  [15]  and  entered  a  new 
phase  for  systematic  coherent  radar  detection  with  Huynen’s  work  on  radar  target 
phenomenology  [16]-[17],  which  led  to  target  decomposition  techniques.  These  techniques 
represent  a  key  theoretical  approach  to  target  recognition  in  current  radar  imagery.  The 
essence  of  target  decomposition  theory  is  to  express  the  average  scattering  mechanism  in 
terms  of  independent  scattering  elements  and  to  associate  each  of  those  with  a  known  physical 
mechanism  or  scatterer  [18].  The  described  association  can  be  on  a  coherent  target  detection 
(CTD)  [17],  [19]-[20],  or  on  a  partially  coherent/incoherent  target  detection  (PCTD)  basis 
[17]-[18].  In  a  CTD  methodology,  the  scattering  matrix  is  fully  polarimetric.  In  a  PCTD 
approach,  however,  the  covariance  or  coherency  matrices  contain  the  polarimetric  information 
of  partially  polarized  scattered  waves.  One  widely  used  CTD  technique  is  Cameron’s 
elemental  decomposition  of  maximum  symmetric  scattering  into  scatterer  classes  of  trihedral, 
dihedral,  dipole,  cylinder,  narrow  diplane  and  quarter  wave  for  each  resolution  cell  [20]. 

There  are  certain  inherent  problems  associated  with  Cameron’s  classification  for  polarimetric 
SAR  (Synthetic  Aperture  Radar)  applications,  including  coarse  segmentation  and  erroneous 
effects  due  to  scattering  segmentation  of  non-coherent  scattering  (under  a  coherent 
assumption).  To  alleviate  these  problematic  effects,  new  methods  have  been  proposed  to 
improve  the  segmentation  process  and  to  make  use  of  coherence  assessment  to  validate  the 
coherent  nature  of  scattering  for  the  target  under  investigation  [21].  The  incoherent  techniques 
(i.e.,  PCTD)  provide  physical  insight  to  the  nature  of  a  partially  coherent  scattered  field  and 
the  degree  of  polarization  through  Kennaugh,  Mueller,  covariance  or  coherency  matrix 
definitions  [17].  These  coherency  indicators  provide  valuable  information  on  the  nature  and 
class  of  the  scatterer  [22]. 
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In  the  present  work,  the  main  objective  is  to  establish  a  relatively  general  and  accurate 
framework  for  modeling  of  EM  propagation  and  scattering  for  airborne  or  spacebome 
polarimetric  applications.  The  theoretical  approach  is  different  from  conventional  techniques 
as  it  combines  a  multidimensional  Fourier  analysis  (i.e.,  decomposition)  of  space  Green’s 
function  and  inverse  analysis  of  scattering  using  the  EM  theory  and  the  principle  of  equivalent 
sources.  Due  to  Fourier  transformation  of  the  exact  source  impulse  response,  coherent 
information  is  included  in  the  formulation.  The  spectral  nature  of  the  analysis  also  provides 
insights  into  the  Doppler  effects  through  the  space-frequency  mapping  mechanism  between 
scatterer  and  observation.  In  addition,  the  forward  nature  of  the  model  is  advantageous  for  the 
purpose  of  target  identification. 

The  other  novel  aspect  of  the  current  method  is  the  classification  of  target  scatterers  through 
the  introduction  of  vector  current  functions.  The  latter,  which  are  based  on  the  concept  of 
scattering-induced  equivalent  sources  in  electromagnetic  theory,  is  close  to  the  roots  of  EM 
scattering  governed  by  Maxwell’s  equations.  To  classify  the  target  scatterer,  the  associated 
vector  current  function  in  the  space  or  spectral  domain  is  decomposed  into  a  set  of  basis 
functions  that  can  be  viewed  as  features.  Various  basis  sets  may  be  used  for  the  required 
decomposition  with  each  set  representing  a  class  of  features  or  sub-features.  Due  to  the 
generic  nature  of  this  formulation,  well-known  analytic  and  numerical  techniques  can  be 
adopted  from  other  areas  (e.g.,  RF,  microwaves)  for  such  characteristic  function  projections 
into  the  desired  feature  space.  The  possibility  of  defining  various  basis  sets  and  number  of 
members  or  basis  functions  for  each  set,  indicates  the  potential  of  controlling  and/or  designing 
the  nature  and  number  of  features  for  target  decomposition.  Such  capability  can  be  used  to 
lower  the  misclassification  probability  associated  with  coarse  classification  or  segmentation 
encountered  through  current  coherent  target  decomposition  techniques.  Furthermore,  study  of 
the  induced  equivalent  current  source’s  behavior  on  a  scatterer,  can  yield  a  more  tangible 
understanding  of  the  scatterer’s  physics.  Finally,  the  source  vector  functions  in  the  present 
formulation  may  be  readily  applied  to  compute  the  associated  EM  field  quantities  (such  as 
coherent  radar  cross  section)  in  both  near  and  far  field  zones. 

An  important  issue  considered  throughout  this  work  and  future  related  works  is  to  identify  and 
examine  various  aspects  of  target  recognition  and  imaging  via  their  contributions  in  the 
analytic  model.  Upon  identification  and  classification  of  these  analytic  entities  or  blocks, 
comparisons  can  be  made  with  the  experimental  results  and  other  techniques  (e.g., 
simulations)  for  the  respective  observables  to  validate  the  model.  Furthermore,  potential 
existing  analytical  tools  suitable  for  characterization  of  a  specific  block  may  be  applied  or 
modified.  One  important  category  is  the  target  scattering  information  with  key  observables 
such  as  the  scatterer’s  EM  signature  and  its  radar  cross  section  (RCS). 

According  to  the  described  goals  and  structure  of  present  methodology,  one  may  outline  the 
overall  objectives  as  follows: 


I.  To  develop  a  general  Fourier-based  formulation  of  electromagnetic  propagation  and 
scattering  for  polarimetric  SAR  applications; 

II.  To  extract  target  scattering  characteristics  from  the  model  and  connect  them  to 
conventional  features,  e.g.,  the  scattering  matrix  and  the  scatterer’s  RCS; 
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III.  To  investigate  inversion  techniques  in  the  spectral  domain  and  potentials  of 
tailoring/optimizing  them  for  certain  types  of  applications; 

IV.  To  study  the  feasibility  of  a  systematic  target  classification  and/or  decomposition  in 
the  spectral  domain; 

V.  To  explore  frequency  domain  selective  techniques  for  adaptive  and  targeted  scattering 
analysis; 

VI.  To  study  the  effects,  significance  and  feasibility  of  an  optimum  collection  path  design 
by  using  the  EM  model;  and 

VII.  To  explore  the  possibility  of  a  multi-incident-angle  analysis  for  target  recognition 
enhancement. 

Achievement  of  these  objectives  is  a  substantial  task.  However,  these  goals  signify  the 
research  units  or  blocks  to  be  investigated  throughout  the  present  and  future  work. 
Formulation  of  EM  scattering  in  the  Fourier  domain  (I),  which  is  at  the  core  of  the 
investigation,  is  presented,  as  well  as  certain  considerations  on  scattering  characteristics  (II) 
and  potential  inversion  techniques  (III).  A  more  in-depth  EM  analysis  of  coherent  scattering 
(e.g.,  RCS)  and  optimization  of  inversion  techniques  will  be  included  in  a  later  report.  Target 
classification  and  decomposition  in  the  spectral  domain  (IV),  potential  adaptive  techniques  in 
the  spectral  domain  based  on  the  state  matrix  and  an  analytic  interpolation  formulation  [23]- 
[24]  to  enhance  target  recognition  (V),  collection  path  design  considerations  (VI)  and,  multi¬ 
look  analysis  of  scatterer  (VII)  using  the  spectral  scattering  model  are  under  study  and  will  be 
reported  in  future  communications. 

Governing  electromagnetic  field  equations  which  are  at  the  core  of  any  EM  field  analysis  are 
described  in  the  next  section.  The  polarimetric  model  is  formulated  in  Section  3,  upon  the 
introduction  of  equivalent  sources  and  derivation  of  the  scattered  EM  field  due  to  the  induced 
sources.  The  general  relation  between  the  described  scattering-induced  sources  and  scattering 
matrix  parameters  are  also  established  in  Section  3.  Section  4  is  devoted  to  the  spectral 
analysis  of  scattering  through  Fourier  decomposition  of  exact  point  source  Green’s  function. 

In  particular,  the  important  system  kernels  are  derived  that  relate  the  polarimetric  SAR 
observations  to  the  spectral  distribution  of  scatterer.  Bandwidth  properties  are  investigated  in 
Section  5.  In  this  section,  a  procedure  to  map  the  windowing  effect  from  space  to  frequency 
domain  is  analytically  detailed  that  provides  an  invaluable  tool  to  investigate  effective 
bandwidth  and  resolution.  It  is  also  highlighted  how  the  target  distribution  in  Fourier  space 
may  be  explored  to  examine  various  radar  incident  and  aspect  angles  by  virtue  of  the 
described  mapping.  The  RCS  derivations  based  on  the  knowledge  of  target  vector  distribution 
are  presented  in  Section  6.  Section  7  includes  concluding  remarks  and  direction  for  future 
research. 
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2.  General  electromagnetic  formulation 


Electromagnetic  fields  in  a  linear,  isotropic  and  inhomogeneous  medium  are  given  by  the 
Maxwell’s  equations: 


Vx 


'VxE' 
,  Mr  , 


-k2£rE  =  -jco/u03xrc 


(1) 


Vx 


^VxH^ 


-k2nr  H  =  Vx 


( J  ^ 

src 

\er  J 


(2) 


In  ( 1  )-(2),  E,H  are  the  complex  electric  and  magnetic  fields  with  time  dependence  eja>l ,  co 
is  the  temporal  frequency,  er  and  /Jr  ( £0  and  ju0  )  are  the  relative  (free  space)  permittivity 

and  permeability,  k  =  is  the  free  space  wavenumber,  c  denotes  the  speed  of  light,  and 

J xrc  represents  all  the  current  sources  within  the  medium  of  interest.  The  differential  operator 
denoted  by  V  is  the  Nabla  or  Hamilton  operator. 

For  scattering  analysis,  it  is  convenient  to  convert  the  above  equations  into  relations  for  an 
equivalent  scattering  problem  in  the  homogeneous  medium  through  introduction  of  equivalent 
induced  polarization  sources  [25]-[26].  Using  the  principle  of  volumetric  equivalence,  (l)-(2) 
can  be  written  as: 


Vx(VxE)-&2E  =  -/<y//0J-VxK  (3) 

V  x  (V  x  H)  -  A:2H  =  V  x  J  -  j(O£0K  (4) 


with  the  equivalent  electric  and  magnetic  current  sources  given  by: 


J  =  J,rr+y^o(fr  “  1)E 


(5) 
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K  =  yo>//0(>, -1)H 


(6) 


Equations  (3)-(6)  are  general  and  provide  a  solid  theoretical  point  of  reference  for  EM 
propagation  and  scattering  modeling  for  radar  (polarimetric  in  particular)  applications.  A 
comprehensive  description  of  scattering  can  be  obtained  by  satisfying  the  boundary  conditions 
and/or  identifying  the  sources  for  each  scenario. 
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3.  Polarimetric  radar  EM  model 


3.1  Field  solution 


To  develop  a  scattering  model  for  polarimetric  radar  applications,  an  electromagnetic  theory 
is  employed.  Maxwell’s  equations  (3)-(4)  are  governing  equations  for  the  total  electric  and 
magnetic  fields  which  include  both  the  incident  (E'  ,H'  )and  scattered 

(E ' ,  H ' )  fields.  For  polarimetric  radars,  a  convenient  scattering  equation  can  be  written  for 
the  electric  field  as: 


V x  (V  x  Ev)  -  k2 E'  =  -jcoii0 J,  -  V x  K, 


(7) 


In  derivation  of  (7),  we  note  that  the  incident  fields  satisfy  the  vector  Helmholtz  equations: 


V2E'  +  A:2E'  =  0 


(8) 


V2H'  +  A:2H'  =0 


(9) 


Note  that  vectors  J  v  and  K  v  in  (7)  are  electric  and  magnetic  currents  that  represent 
equivalent  induced  sources  due  to  scattering. 

Although  modeling  of  the  above  equivalent  sources  introduces  a  challenge  for  certain 
applications,  the  approach  provides  a  robust  theoretical  tool  for  the  analysis  of  scattering. 
Another  attractive  feature  is  the  possibility  of  incorporating  medium  effects  in  definition  of 
the  equivalent  sources  (i.e.,  in  equations  (5)-(6)). 

A  general  solution  to  Maxwell’s  equations  (e.g.,  (7))  for  the  scattered  field  at  the  observation 
point  r ,  can  be  found  by  taking  advantage  of  the  well-known  three-dimensional  Green’s 
function, 
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(10) 


G(r)  = 


4^r|r| 


as: 


E«'(r,®)  =  —(*’  +  VV)[G  •  J,]-  V  x  [G  •  K ,]  (1 1) 

JCOS  o 


where  the  asterisk  (*)  denotes  three-dimensional  space  convolution  and,  as  before,  the  time 
dependence  is  eJ0>' . 

Making  use  of  the  properties  of  the  free-space  Green’s  function,  the  far-field  expression  (i.e., 
the  radiation  gauge)  becomes: 


E«'(r,ffl)=-4—  (G*  J,)r -(Vx  (G.K,)  ),.  (12) 

jas  o 


where  the  subscript  T  denotes  the  transversal  component  in  a  plane  orthogonal  to  the 
direction  of  radiation  or  line  of  sight  (LOS).  The  radiation  field  in  (12)  can  also  be  written  as: 


Erec'vr  (r,co)  =  jk  {  7[(G  *  J  v )  x  r ]  x  r  +  (G  *  K  v )  x  r }  (13) 


where  rj 


=  k/  =  \Mo 

/0)Sn 


is  the  intrinsic  impedance  of  the  medium  (free  space)  and  r  is 


the  unit  vector  in  the  direction  of  radiation  (i.e.,  observation). 
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3.2  Equivalent  sources 


There  is  no  constraint  on  the  type  of  induced  current  sources  in  (13);  in  general  they  can  be 
volumetric.  For  scattering  analysis  in  radar  applications,  however,  it  is  convenient  to 
formulate  the  problem  in  terms  of  surface  distributed  sources.  Valid  electromagnetic 
formulations  of  such  boundary-value  problems  provide  a  comprehensive  picture  of  scattering 
[25]-[28].  The  main  objective  of  these  formalisms  is  to  derive  the  overall  field  (initial  and 
induced  or  scattered)  by  solving  respective  integro-differential  equations  for  fields  and  source 
distributions  when  the  required  boundary  conditions  are  enforced.  It  is  crucial  to  use  proper 
boundary  conditions  for  integro-differential  formulations  in  order  to  produce  unique  field 
solutions.  In  particular,  care  should  be  taken  in  defining  equivalent  sources  when  fields  are 
non-vanishing  infinitesimally  close  to  the  boundary  surface.  The  latter  may  cause  erroneous 
results  even  for  infinitesimally  thin  layers  of  p.e.c.  (perfect  electric  conductor)  or  p.m.c. 
(perfect  magnetic  conductor). 

As  was  addressed,  the  fundamental  theorem  of  electromagnetics  generally  known  as  the 
equivalence  principle  is  a  powerful  concept  that  is  rich  in  physical  insight  and  provides 
theoretical  ground  for  relating  the  fields  to  scattering  through  definition  of  equivalent  induced 
sources.  The  concept  of  equivalent  or  substitute  sources  is  a  well  known  approach  that  can  be 
traced  back  to  Huyghens  [29],  This  concept  was  further  developed  by  Schelkunoff  [30], 
whose  equivalence  principle  asserts  that  use  of  surface  equivalent  sources  derived  from  the 
overall  field  values  produces  true  fields  at  all  points  outside  the  surface  and  null  fields  inside. 
The  Stratton-Chu  formulation  [3 1]  provides  the  same  result  when  it  is  used  to  express  the 
electric  and  magnetic  fields  in  terms  of  surface  integrals  [32],  One  can  apply  the  described 
surface  equivalence  principle  for  target  scattering  modeling  when  the  necessary  boundary 
conditions  are  satisfied. 

Assume  S  is  the  scattering  surface  of  the  target  or  scatterer  with  electromagnetic  fields  E,0, 
and  H/0,  uniquely  defined  at  this  surface.  Employing  Schelkunoff  s  equivalence  principle,  we 
may  write  the  electric  and  magnetic  currents  as: 


J. v  (r.v ) =  n  v  (rv ) x  H„„  (rv ) 

(14) 

K  ,v  (r, )  =  -  n,  (r, )  x  E,„,  (rv ) 

(15) 

where  nv(rv)  denotes  the  unit  (outward)  normal  vector  to  the  surface  S  at  scattering  point 

rv .  The  total  fields  E,0,  and  Hlnl  at  any  point  of  observation  space  can  be  considered  as  the 

sum  of  the  incident  fields,  i.e.,  the  electromagnetic  fields  when  there  is  no  scattering,  plus 
secondary  or  scattering  fields: 
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E„„(r)  =  E'(r)  +  E'(r)  (16) 

H„(r)-H'(r)  +  H'(r)  (17) 

In  general,  one  can  set: 

E/0/(r)  =[i  +  S£(r)]-E'(r)  (18) 

H,„,(r)  =  [I  +  Sw(r)]-H'(r)  (19) 


where  SE  and  SH  are  dyadic  quantities  representing  general  electric  and  magnetic  field 
scattering,  and  I  is  the  unitary  dyadic. 

In  theory,  one  may  define  the  equivalent  sources  at  any  arbitrary  surface  S  in  space  as: 


J, (rt)  =  ns(r,)x{[I  +  Sw(rv)]-H'(r,)}  (20) 

K ,  (r, )  =  -  h ,  (r, )  x  { [I  +  S  ,  (r, )]  •  E'  (r, ) }  (21) 


The  sources  in  (20)-(21)  produce  true  fields  outside  the  surface  S  and  a  null  field  inside. 
Defining  the  general  scattering  from  electric  incident  to  magnetic  scattering  YHE  with: 

H1  =  Yhe  -E'  (22) 
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one  may  recast  (20)  as: 


J,(r,)  =  n,(r,)x U  x E'(r,)  +  Y„/;(r,).E'(rv) 


(23) 


In  (23),  r(  is  a  unit  vector  specifying  the  propagation  direction  of  the  incident  EM  wave 
(E',H') .  Application  of  (21)  and  (23)  in  (13)  yields  an  expression  for  the  scattered  or  the 
echoed  field  at  the  space-time  observation  point  PKC(r,t)  . 

At  this  stage,  no  particular  form  is  assumed  for  S;.. ,  Sw  or  Y  H,  which  are  simply  reflecting 
the  general  scattering  mechanism,  i.e.,  the  difference  between  the  overall  field  and  the  initial 
or  incident  field.  As  will  be  discussed,  however,  these  quantities  can  be  directly  related  to  the 
polarimetric  scattering  matrix  for  a  backscattering  scenario. 


3.3  Radar  field  model 


By  defining  the  right-handed  radar  system  unit  vectors  (h  ,v,  r, ) ,  the  known  incident  field 
polarization  at  any  point  in  space  is  represented  by: 


«/„■  (r)  =  cos(^(r))h  +  sin(y/(r))v 


(24) 


where  y/(r)  denotes  the  polarization  angle.  That  is: 


E'(r)  =  £'(r)u,,v(r) 


(25) 


Use  of  (25)  in  (21)  and  (23)  yields: 
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Kv  =-n4.  x[(i  +  S£)-uAl,]i?' 


(26) 


J  v  =  n,  x 


-L*K  +  X 
v 


HE  UAv 


E‘ 


(27) 


All  quantities  in  (26)-(27)  are  computed  at  the  scattering  point  (rv)  .  The  vector  terms  in  (26)- 
(27): 


W  =  -n,x[(f +  §,)•«*,] 


(28) 


F  =  n,  x 


r, 

—  xu 
V 


hr 


+  Y, 


HE 


u, 


(29) 


can  be  viewed  as  (and  hereafter  referred  to  as)  a  polarimetric  radar  target  function. 

Some  considerations  must  be  made  before  completing  the  present  discussion.  The  target 
functions  in  (28)-(29)  represent  the  normalized  equivalent  sources  due  to  scattering.  As  a 
result,  these  vectors  contain  the  coherent  scattering  information  and  can  be  analyzed  to 
estimate  the  scattering  matrix.  In  fact,  for  a  backscattering  scenario  where  targets  are  modeled 
as  reflectors,  the  overall  scattering  mechanisms  in  (28)-(29)  reduce  to  conventional 
polarimetric  scattering  quantities.  Since  all  components  of  these  equations  are  known,  except 
the  scattering  quantities,  target  reflection  functions  (or  equivalent  source  functions)  can  be 
connected  to  the  scattering  parameters.  The  latter  provides  a  means  for  studying  scattering 
characteristics  and  other  target  properties,  such  as  the  coherent  RCS,  within  the  same  analysis. 

Although  possible  to  handle  analytically,  throughout  the  rest  of  this  analysis  magnetic  current 
sources  (28)  will  be  neglected  (i.e.,  we  take  the  electric  reflector  assumption).  It  will  also  be 
useful  to  consider  the  following.  If  the  scattering  matrix  for  the  electric  field  in  a  transversal 
plane  (or  slant  plane)  is  defined  by: 
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(30) 


(31) 


’hh  °/»v 


V^vA 


E ' 

vc»  y 


=  S  •  E' 


and  if. 


(Y, 


hh  *  hv 


1  vv  J 


\{  p  i  \ 

Lh 


=  Y-E' 


one  gets  (see  (8)-(9)), 


Y  =  -  <J2  -S 
V 


(32) 


where 


0  f 
1  0, 


(33) 


is  a  Pauli  spin  matrix. 


Making  use  of  (13),  (27)  and  (29)  gives  the  expression  for  the  scattered  field: 


Erefvr(r ,co)  =  jkr]E'[(G  *F)xr]xr 


(34) 
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The  convolution  in  (34)  is  performed  over  the  coordinates  of  the  scatterer.  Hence,  the  integral 
form  of  (34)  using  (10)  can  be  written  as: 


Erec. (r,  (o)  =  jk  77  JJj£'  (F  x  r  x  f ) 


An  |  r  -  r,  | 


drr 


(35) 


The  integration  in  (35)  is  a  three-dimensional  volume  integral  whereas  the  equivalent  sources 
are  two-dimensional  surface  currents.  This,  however,  does  not  introduce  inconsistency  or  loss 
of  generality  if  the  currents  are  properly  defined.  For  instance,  if  the  current  surface  is  defined 
by  a  normal  vector  n(r,)  and  coordinates  (x„  ,xl2,xn),  the  3-dimensional  source 

distribution  becomes  J(xn,xl2,xn)S(xn  -xm).  In  the  latter,  £(x)  represents  the  Dirac 
delta  function. 

For  convenience,  henceforth,  the  target  or  scatterer  coordinates  are  denoted  by  r ,  the 
observation  coordinates  by  u ,  and  the  unit  radial  vector  by  u .  Moreover,  the  transversal 
normalized  current  Fxrxr  (or  F  x  u  x  u  )  is  represented  by  Fr .  Accordingly: 


Ere/r(u  ,co) 


e-jk\n-r\ 

|u-r| 


dr 


(36) 
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4.  Spectral  domain  modeling 


Multidimensional  Fourier  analysis  of  the  backscattered  field  (36),  provides  valuable  insight  to 
the  nature  of  the  scatterer.  Prior  to  such  a  transformation,  certain  considerations  are  made. 

In  the  present  formulation,  the  centre  of  spatial  coordinates  is  assumed  to  be  the  centre  of 
target  or  scatterer  coordinate  system,  i.e.,  the  target  is  centred  at  the  origin  of  the  coordinate 

system  r  =  (x,  y,  z)  =  (0,0,0) .  The  incident  field  ( E ' )  effects  are  evaluated  at  the  target 
position  r  whereas  the  signal  or  wave  originates  from  the  radar  position  u .  In  order  to 
account  for  this  delay,  one  may  consider  the  incident  field  (at  the  target)  as: 


E'  =  E\o))e-jk'r-ul 


(37) 


where  E'(a> )  is  the  original  transmitted  waveform’s  spectrum.  Thus, 


JJJ  F,(r  ,a>)e-^--dr 


~j2A|u-r| 


An 


u  -r 


(38) 


Furthermore,  assuming  the  amplitude  variations  associated  with  Green’s  function  (far  field) 
for  a  target  area  are  negligible,  one  may  use: 


Ere/ar(u,«y)  =  Jjf  F,(r,fi>)  e~J2k^  dr 


(«  =1  u  |)  (39) 


4.1  Fourier  analysis 


The  general  three-dimensional  (observation)  space  Fourier  transform  of  (36)  is  given  by: 


14 


DRDC  Ottawa  TR  2004-1 15 


where  k„  =  (km,kuy,ku:)  is  the  spatial  frequency  or  wave  vector,  and 


-j2kju]+ul+u] 

g  e~Jku*"x  e~jKyUy  <TAa  dux  duv  du2  (41) 

V  2  2  2  ' 

ux  +  uy  +uz 

Fr(k,o>)=  JJJ  Fy  (r,&>)  e~Jk*x  e~jkyy  e~'k‘2  dxdydz  (42) 

The  integration  in  (41)  is  performed  over  the  observation  space  whereas  the  integrations  in 
(38)  and  (42)  are  performed  over  the  target  or  scatterer  coordinates.  As  can  be  seen  from  (40)- 
(42),  the  three-dimensional  space  frequency  mapping  is  represented  through  kw  =  k  . 

In  order  to  relate  the  observations  (in  either  the  spatial  or  the  spectral  domain)  to  the  Fourier 
characteristics  of  the  scattering  source,  a  proper  Fourier  decomposition  of  the  Green’s 
function  in  (38)  is  required.  It  is  possible  to  handle  the  required  derivations  by  using  the 
principle  of  stationary  phase  (PSP),  as  described  in  Appendix  B.  To  explore  more  physical 
insights  and  a  potential  slant  plain  analysis  for  SAR  applications,  however,  the  exact  Fourier 
analysis  is  carried  out  on  the  point  source  Green’s  function  in  separate  steps  for  the 
transversal  plane  and  the  longitudinal  (or  vertical)  dimension. 

The  source  response  function  in  (38): 


An 


G„(  r)  = 


47r|r| 


satisfies  the  scalar  Helmholtz  equation: 


(43) 


(V2+Ak2)Gu(r)  =  -S(r) 


(44) 


Hence,  the  three-dimensional  Fourier  transform  satisfies  the  following  equation: 
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(*,’  +  *,’  +k‘-4k1)G,(k„kr,k,)  =  l 


(45) 


The  polar  spectrum  of  the  Green’s  function  can  be  obtained  by  the  inverse  Fourier  integration 
over  the  complex  k:  plane.  By  applying  the  residue  theorem,  one  gets: 


Gup(kx,ky,z) 


jtxp(j^4k2  -  p2  \z\t 

2yfik 2  -p2 


where  ( p,(p )  are  the  spectral  polar  coordinates  of  (kx,ky),  i.e., 


(46) 


kx  =  pcos{(p) 

(Ala) 

ky  =  psin(^) 

{Alb) 

kx  +  ky  =  p2 

(47c) 

Next,  we  apply  (46)  in  (38)  when  making  use  of  the  generalized  Parseval’s  theorem  to  obtain: 


S*(» ,*)  =  E\co)  J/J  Fr(kx,ky,z,co) 


dkdkdz 

■*  y 


(48) 


where  the  normalized  received  signal  is  defined  by: 


Srec(u,<y)  =  ~~T~  E* 

krj 


r(u  ,co) 


(49) 
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Hence  for  uz  -  z  >  0,  which  introduces  no  loss  of  generality, 


S„,c(u,a>)  =  £>)JJ  F ,,{kx,ky,kz,co) 


exp(;V4A:2  -/?2  -Jktu,  -jkyUy 


dkxdky 


(50) 


with 


k:  =  ^  4k2  -  p2 


(51) 


Consider  a  general  observation  point,  or  in  a  SAR  scenario,  the  synthetic  aperture  coordinates: 


u  =  ( ux,uy,uz )  =  (r„  cos (</>),ru  sin((Zi),wz) 


(52) 


where  ru  and  (f>  are  polar  coordinates  of  observation  or  synthetic  aperture.  Making  use  of 
(52)  results  in: 


Srec(u,<y)  =  E'(o))  JJ  F r{kx,ky,kz,co) 


expl 

hk2-p 2  mJ 

V 

4 k2  -p2 

exp {-jp  ru  cos(</>  -cp))pdpd<p 


(53) 


where  kx,ky,kz  are  defined  by  (47)  and  (51). 


One  can  recast  (53)  as: 
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~  (“>  ®)  =  (*>)  — W4^2  4  —  [  *V  (*, .  ■ *„ ,  , 


fik: 


,co)*exp(-jpru  cos(^))](f>) 

(54) 


In  (54)  the  asterisk  symbol  (*)  denotes  the  one-dimensional  angular  convolution  ( ^ 
dimension)  and  kx,ky,k.  are  defined  by  (47)  where  (p  is  replaced  with  <f> . 

Equation  (54)  relates  the  spectral  distribution  of  the  polarimetric  target  function  Fy.  to  any 
arbitrary  observation  point  (or  synthetic  aperture)  u  in  space. 

For  moving  platforms,  as  was  addressed  earlier,  space  frequency  mapping  to  the  observation 
spectral  domain  takes  place  in  (54)  (also  see  (40)).  More  specifically,  taking  the  Fourier 
transform  of  both  sides  in  (54)  with  respect  to  u, : 


^reX^uz^x^y^)  &  (®)  J H====  =  ^(Kz  ~  p1  )[  Fy •(*,  Xy  ,  k.  ,  CO)  *  QXp(~  j p ^  COS(0))] 

(55) 


Taking  advantage  of  the  analytic  properties  of  the  delta  function: 


\g(p)d(f(p)-a)dp 


gjp) 
df  I  dp 


P=P» 


(where  f(p0)-a  =  0) 


(56) 


we  find: 


XK:>r„  cos(^),rH  s\n(<j>),co)  =  EX(o)[%(kx,kyX^co)*^p(-jkprucos(<f>))}^  (57) 


where  now, 
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(58) 


K  =  Kz 


K  =  ^4k2  -  kj  cos(^)  (59a) 

ky  =  ^4k2-kJ  sin(^)  (596) 

with 


*,=>2-C  <6°) 


Equations  (58)-(60)  indicate  that  the  mapping  from  the  observation  spectral  domain  to  target 
spectral  domain  is,  in  general,  non-linear.  The  space  frequencies  in  (58)-(59)  are  located  on 
the  surface  of  a  sphere  of  radius  2k  in  the  spectral  domain,  known  as  the  Ewald  sphere. 

Considering  the  Fourier  expansion  of  the  nonlinear  phase  term: 


exp(-/scos(^))  =  ^Hj2)(s)eja(^’r/2)  (61) 

a 


where  Ha(1)  (s)  is  the  Hankel  function  of  the  second  kind  and  order  a  ,  one  may  transform 
(57)  to  the  spectral  domain  (a)  of  angular  dimension  <f>  (which  can  also  be  considered  as  an 
aspect  angle)  as: 


S rec(ku2,cc,ru,CQ)  =  E'(co)  F,1 


{<*) 


ikuz,a,co)Ha(1\kr)e 


-Jan  j2 


(62) 


The  target  function  Fourier  transform  is  defined  via: 
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(63) 


¥r(a\kuz,cc,co)=  Jf T(kx,ky,ks,co)e-Ja*  d<j> 
*» 


Fy.  {kx  ,ky,kz,co)  =  Yj  ^ra)  (Kz  >a>  <°)  e'a* 


(64) 


with  (kx,ky,k.)  defined  by  (5 8)-(59)  and,  </>0  representing  the  angular  synthetic  aperture 

domain.  It  is  within  this  domain  that  orthogonality  of  Fourier  components  e~ja *  holds. 
Accordingly,  for  <j)  e  ( -n,n )  ,  a  takes  on  integer  values. 

Equation  (62)  (or  (57))  provides  an  analytic  means  to  recover  the  target  information  contained 
in  F7  from  the  observations  made  through  STO.  A  variety  of  methods  may  be  used  to 

perform  the  required  inverse  processing  which  are  optimal  and  suited  for  certain  applications. 
A  relatively  simple  and  general  matched  filtering  technique  to  implement  the  deconvolution 
needed  for  inverse  processing  will  be  addressed. 


4.2  2-Dimensional  slant  plane  Fourier  analysis 


As  a  typical  SAR  scenario,  consider  the  imaging  of  the  ground  plane  by  a  radar-carrying 
vehicle  flying  at  a  constant  height  z  =  Zh.  Although  altitude  variations  from  the  ground  plane 

z  =  0  are  always  present  in  practical  operations  (as  included  in  the  model  developed),  those 
effects  may  not  be  of  interest  for  certain  applications.  An  appropriate  model  for  such 
applications  can  be  derived  from  the  three-dimensional  formulation  by  imposing  the 
associated  conditions.  We  perform  the  required  derivations  as  follows. 

The  target  distribution  at  the  ground  level  z  =  0  (i.e.,  ignoring  the  terrain  variations)  can  be 
written  as: 


Fy  (*,  y,  Z,  co)  =  F,  (x,  y,  co)S  (z) 


(65) 


Applying  the  three-dimensional  Fourier  transform  of  (65)  in  (57),  one  gets: 
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Snc (kU2 , ru  cos {</>), ru  sin (</>), co)  =  Ei (co) [  fr (kx ,ky,co)* exp (~jkp  ru  cos(</>)) ]  w  (66) 


The  inverse  Fourier  transform  of  (66)  with  respect  to  uz  (kuz )  yields: 


Snc(uz,ru  cos (<f>),ru  sin (</>),co)  =  E' (co) 


\dkuz  eJk”"!  [ Fr ( kx,ky,a >)  * exp (~jkp  ru  cos(^)) ] w) 


(67) 


where  k  p  is  given  by  (60).  Hence, 


(68) 


Making  use  of  (68)  for  the  observation  set  at  uz  =Zh,  results  in: 


Ic  {jIIc  I  2  2  ^ 

Sret (Zh , ru  cos (<j>), ru  sin(^), co)  =  E‘ (co)  J  .  P  p==e'j^k  ~kp  z‘  [ F> (kx ,ky,co )  * exp (~jkp  r„  cos(<j>)) ] w 

V4* 


(69) 


with, 


kx  =  cos((Z)) 


^  =  kp  sin(^) 


(70) 


(71) 


Furthermore,  the  angular  Fourier  transform  described  in  (63)-(64)  may  be  used  to  transform 
(69)  into: 


DRDC  Ottawa  TR  2004-1 15 


21 


(72) 


Snc(Z„,a,ru 


“7 yj^k^-k 

co)  =  E\co )  \kpdkp  ~  =  ¥;(a) (kp,a,co)  Ha{1\kpru )  e'^2 

V4*  -kp 


As  can  be  seen  from  equation  (72),  extraction  of  the  ground  target  information  ¥r  from  the 

observation  Sn,c  \z=zh  's  more  complicated  than  in  (62)  and,  in  principle,  involves  solving  the 

integral  equation.  One  should  note,  however,  that  the  simplicity  offered  in  (62)  is  due  to  a  full 
three-dimensional  spectral  domain  treatment  of  scattering,  i.e.,  an  extra  layer  of  Fourier 
transformation  in  the  z  direction.  In  theory,  the  system  defined  by  equation  (72)  can  be 
viewed  as  a  shift-varying  system  [33]  with  a  known  analytical  kernel: 


Y{co,co  )  =  E\co)  =■  (73) 

V4*-V 


where  the  ground  plane  equivalent  temporal  frequency  m  is  defined  via: 


(74) 


Setting: 


(75) 


the  system  equation  can  be  written  as: 


Sn!lXZh,a,ru,ct))=  jr(Zh,coK,ct))  S g(a,rt,a>g)da)g 


(76) 
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Technically,  the  ground  plane  signal  (75)  can  be  recovered  by  using  the  inverse  of  the  system 
kernel  function  (73): 


S  ^a,ru,(Og)=\Srec{Zh,a,ru,o})T-\Zh,o)g,o})d(o 


Application  of  (77)  with  (75)  yields  the  desired  target  function.  One  should  note  that  the 
required  deconvolution  to  recover  the  ground  plane  signal  could  also  be  achieved  through 
matched-filtering  of  the  system  kernel  scheme. 
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5.  Fourier  properties  and  inversion  considerations 


The  model  developed  here  can  be  used  to  examine  various  aspects  connecting  the  target 
characteristics  to  observations  from  both  electromagnetics  and  system  points  of  view.  One 
important  application,  however,  is  to  extract  the  target  information  from  observations. 
According  to  equations  (62)  (or  (57))  and  (72)  (or  (62)),  the  latter  involves  inverse  processing 
utilizing  techniques  ranging  from  matched  filtering  methods  (as  will  be  discussed)  to  more 
sophisticated  integral  equation  based  methods  such  as  the  method  of  moments  (MoM). 

To  perform  the  inverse  processing  and  conduct  any  spectral  domain  considerations  (e.g., 
spatial  resolution),  bandwidth  limitations  associated  with  factors  such  as  finite  aperture  and 
target  distribution  must  be  investigated. 


5.1  Bandwidth  properties 


Bandwidth  characteristics  of  the  received  SAR  signal  can  be  studied  by  making  use  of  the 
windowing  properties  of  phase-modulated  (PM)  signals  [34]-[35],  These  properties  are 
essentially  associated  with  the  space-limited  nature  of  the  Green’s  function  connecting  the 
target  (scatterer)  and  observation  spaces.  Since  the  Green’s  function  and  its  spherical  phase 
front  are  space-limited  or  windowed,  the  Fourier  transform  of  the  Green’s  function  is  also 
windowed  in  the  space-frequency  domain.  To  formulate  the  described  windowing  behaviour 
for  bandwidth  estimation,  one  can  impose  the  space  windowing  conditions  when  deriving  the 
instantaneous  frequencies. 

We  will  derive  a  general  relation  to  map  the  respective  spatial  and  spectral  windowing.  Prior 
to  that,  however,  the  following  should  be  noted.  The  Fourier  decomposition  of  the  Green’s 
function  in  the  present  formulation  represents  all  possible  components  including  evanescent 
modes.  For  SAR  applications,  the  observation  is  normally  within  the  far  field  region  of  the 
target  where  the  evanescent  or  transient  modes  are  negligible.  Hence,  considering  the 
characteristic  equation  for  the  poles  of  the  Green’s  function  (Ewald  sphere): 


kx  +  ky  +k22  -4k2  =  0 


(78) 


the  overall  support  band  in  the  frequency  domain  (infinite  aperture  and  target  distribution)  is 
2k  ,  or  equally,  (kx,ky,k.)  e[-2k,2k]. 
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The  generalized  instantaneous  space  frequency  vector  can  be  obtained  by  computing  the 
gradient  (V)  of  the  phase  function  'F(r)  with  respect  to  space  coordinates: 


k-W=W(r) 


(79) 


Equation  (79)  provides  an  effective  tool  to  evaluate  the  spectral  domain  properties  of  domain 
limited  field  signals.  One  can  show  that  (79)  satisfies  the  generalized  stationary  phase 
equation: 


TO-kw(,)-r]  =  0  (80) 

that  is  (in  scalar  format)  used  to  determine  the  phase  center  of  the  phase  function  [35]. 
Application  of  (79)  for  the  spherical  phase  front  associated  with  the  Green’s  function: 


T(u,r)  =  2£|u-r 


(81) 


yields  the  instantaneous  frequency  for  spectral  domain  considerations.  Spectral  domain 
windowing  properties  are  examined  throughout  the  following  analysis. 

Consider  the  SAR  operation  for  ground  plane  target  analysis  detailed  in  Section  4.2.  In  this 
typical  case  of  operation,  the  observation  undertakes  a  general  motion  at  a  certain  altitude 
z  -Zh .  Furthermore,  we  assume  the  target  has  a  limited  distribution  within  a  radius  of  r0 , 

i.e.,  r  =  ^[x^+yT  <  r0  (Figure  1). 


Using  (79)  and  (81),  one  finds  the  ( x,y )  plane  instantaneous  space  frequencies: 
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with  associated  polar  coordinates  ( p ,  cp)  given  by: 


Define  the  centre  slant  angle  associated  with  the  centre  of  target  distribution  (centre  incident 
angle),  as  given  in  Figure  1(a): 


Inspection  of  Figure  1(b)  reveals  that  the  spectral  polar  angle  (p  defined  in  (84)  is  windowed 
due  to  the  finite  target  distribution  as: 

<P^<p  o  (87) 

where, 
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<Po  =  sin”'  ( 


=  nj  2 


r0  ^  r« 

otherwise  (88) 


Making  use  of  (85),  we  find  the  polar  phase  center  for  any  given  observation  point  (ru) : 


Po  =  2k  -  =====  2k$m(6z) 


(89) 


with  a  spread  that  can  be  approximated  by: 


A p  »  2&cos2(i9;)sin(^0) 


(90) 


Thus  the  windowing  function  in  the  spectral  domain  for  a  given  observation  point  ( ru )  is: 


Wr(p,<p,co)  = 


for  |  p  —  2k  sin(6»_, )  |  <  2k  cos2  (0Z ) sin(^0)  and  \(p\<(p0 
otherwise  (91) 


For  a  more  particular  and/or  typical  case  of  operation,  one  may  consider  a  scenario  when  the 
observation  is  performed  along  a  straight  line.  In  this  case,  the  observation  vector  is 
represented  by: 


u  —  {2C,uy,Zh) 


(92) 
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where  X  denotes  the  constant  ground  range  associated  with  the  center  of  target  distribution 

and,  -L<uy  <L  represents  the  synthetic  aperture  dimension.  Accordingly,  L  is  the 

synthetic  aperture  length.  The  instantaneous  radial  frequency  for  the  present  scenario 
becomes: 


P 


2k 

^(X-x)2  +  (uy  -y)2  +Zh2 


(93) 


Define  the  instantaneous  centre  incident  angle  associated  with  any  observation  point  ( r  )  as 


6 ,  =  sm~' 


(94) 


where. 


(95) 


The  spectral  windowing  function  for  any  observation  point  ( ru ),  is  given  by  (89)-(91)  when 
definitions  (94)-(95)  are  applied.  Hence,  moving  along  the  uy  direction  from  -  L  to  +  L , 
i.e.,  uy  e[-Z,Z]  : 


2k sinC^ (uy )) - 2k cos 2 (6, (uy )) sin(^0)  <  p  <2ks\n(G:(uy))+2kcos2{e:(uy))sm((pQ) 

(96) 


In  order  to  determine  the  window  function  associated  with  (96),  consider  the  functions: 
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(97) 


— Q  (y)  =  —[sin O')  -  cos2  (r)sin(^0 )]  =  cos(/)[l  +  2  sin(r)sin(^0 )] 
dy  dy 


— Q  (y)  =  — [sin(^)  +  cos2  ( y )  sin(^0 )]  =  cos(;k)[1  -  2  sin(^)  sin(^0 )]  (98) 

dy  dy 


It  can  be  seen  that  — ^  is  positive  for  0<  y  <njl ,  while  —  becomes  negative 
dy  dy 

at: 


sin(r,nax)  = 


1 

2sin(^0) 


(99) 


The  maximum  function  value  associated  with  ymm  is: 


n+(r«) 


l  +  4sin2(ff0) 
4sin(^0) 


(100) 


Equation  (97)-(100)  indicate  that  within  the  interval  [y0,yL]  (D'o>//J  e  [0,tt/2] ),  one  can 
write: 


sin(/0 )  -  cos2  (y0 ) sin(^0 )  <  Q_(y)  <  sin(/,)-cos2(^)sin(^0)  (101) 
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sin(/o )  +  cos2  (y0 ) sin(^0 )  <  Q+  (y)  <  sin (/, )  +  cos2  ( y , ) sin(^?0 ) 

sin(^o)  +  cos2(70)sin(^0)  <  Q+(y)  <  —  4sin 

4sin(^0) 


>for  Y,  <  rmax 

J°r  Y  l  >  Ym  ax 
(102) 


Taking  advantage  of  ( 1 0 1 )-( 1 02),  the  support  band  associated  with  (96)  is  as  follows: 


2k sin(0.o )-2k cos 2 (0lO ) sin(^0 )  <  p  <  2ksm(0zl)  +  2kcos2(0;,)sm((po )  ,far  0.,  <ymn 

J°r  6-j.  >Ymsx 

(103) 


2k  sin(#.0  )-2k  cos2  (0zO )  sin(^0 )  <  p  <  l+4sin 

4sin(^0) 


The  respective  windowing  function  becomes: 


Wr(p,<p,a))  = 


1 

0. 


for  p  satisfying  (103)  and\<p\<(pQ 
otherwise 


(104) 


where  ymm  is  defined  by  (99)  and: 


0:Q  =  tan_,(-^-) 

Z h 


(105) 


0zl  =  tan  1 


Jx2+L2 


(106) 
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One  last  note  in  the  context  of  windowing  characteristics  of  phase-modulated  signals  is  the 
windowing  properties  when  the  angular  Fourier  transform  ((63)-(64))  is  used.  In  this  case,  the 
support  of  the  window  for  the  angular  direction  (j>  in  the  respective  spectral  domain  a  must 
be  determined.  The  latter  can  be  done  by  applying  the  instantaneous  frequency  relation  (79)  in 
the  angular  direction  and  the  phase  definition: 


exp (~jpru  cos {$))  =  exp(-y'T) 


(107a) 


as, 


a 


inst 


d(j> 


-- Pru  sin(^) 


(1076) 


Thus, 


\a\<prusm{cpQ) 


(108) 


As  was  asserted  earlier,  derivations  in  this  subsection  can  be  effectively  used  to  perform 
spectral  analysis  of  the  polarimetric  radar  waveform  for  different  purposes.  One  typical 
application  (as  detailed)  is  the  effective  bandwidth  estimation.  Another  potentially  interesting 
application  is  the  exploration  of  acquired  data  along  the  paths  corresponding  to  different 
incident  angles. 

Consider  the  description  of  instantaneous  frequencies  given  by  (82)-(83),  (85)  and  a  similar 
description  for  kz  that  can  be  easily  derived  in  the  same  fashion.  Making  use  of  the  above  and 
(86),  one  can  derive  a  similar  relation  for  the  incident  angle  in  spectral  domain,  i.e., 
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Various  paths  (corresponding  to  certain  incident  and  aspect  angles)  within  the  3-dimensional 
spectral  distribution  of  the  target  scatterer  may  be  explored  in  the  spectral  space 
(kx ,  ky ,  kz )  by  using  ( 1 09)  and  (84). 


5.2  Inversion 


As  was  addressed  earlier,  recovery  of  the  target  characteristic  vector  function  requires 
inverse  mapping  according  to  the  derived  scattering  model  (e.g.,  (57),  (69)).  Such  inverse 
processing  normally  involves  deconvolution  of  the  known  system  kernel  or  reference  function 
and  may  be  performed  through  a  variety  of  analytic  and/or  numerical  methods.  One 
potentially  strong  candidate  is  the  method  of  moments  (MoM)  or  its  variational  form, 

Galerkin 's  Method  [36]-[39].  Moment  methods  are  suitable  and  numerically  efficient  for 
computations  involving  known  analytical  kernels  [36],  We  do  not  intend  to  discuss  the 
analytical  details  of  such  formulations.  In  general,  the  procedure  is  to  solve  the 
inhomogeneous  system: 


L  (f)  =  g 


(110) 


where  L  is  the  system  operator  which  can  be  viewed  as  the  system  kernel,  g  is  the  known 
quantity  (e.g.,  observation),  and  /denotes  the  unknown  function,  e.g.,  scattering 
characteristics  to  be  determined.  Expanding  /  into  a  set  of  appropriate  basis  functions: 


/=£«*/. 

rn 


(111) 


and  defining  a  proper  projection  space  {/?,.}  and  moments, 


M =<  w/H,v,.> 


(112) 


gives: 
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(113) 


2>.  <L/™  ’A  >=  > 

m  m 


The  resultant  linear  system  is  solved  for  the  unknown  weightings  am .  In  (1 13),  the  system 
operator  L  is  assumed  to  be  self-adjoint.  Proper  selection  of  the  basis  set  {fm  } ,  projection 
space  {pj} and  mechanism  (  u,v  )  (generalized  inner  product)  is  important  to  yield 
convergence  and  accurate  results.  Choosing  the  projection  space  {/?,  }  to  be  the  same  as  the 
basis  functions  space  {fm} ,  results  in  the  variational  Galerkin’s  method. 

A  formalism  based  on  the  method  of  moments  can  be  developed  and  tailored  for  the  present 
SAR  inverse  processing  to  recover  the  polarimetric  target  characteristics  through  knowledge 
of  the  analytic  system  kernel.  One  can  explore  various  basis  vector  space  design  (under 
investigation  and  to  be  communicated)  representing  scattering  decomposition  classes  for  the 
purpose  of  target  identification. 

Relations  (62)  and  (75)  suggest  that  the  target  function  in  the  spectral  domain  may  be 
recovered  by  simply  dividing  the  observation  spectrum  by  the  system  kernel  or  reference 
function: 


F,(a>  (km  ,a,a>)  =  SK/~'  (kp,a,  ru ,  a)  Srec.  (k,a ,  a,  ru ,  to) 
F/o)  (kp ,  a)  =  Sgre (kp  ,a,ru )  Sg  (a,  ru ,  cog ) 


where. 


^  „-W2. 


SKf  (kp,a,ru,co)  =  E'(a>)  Ha  (kpru )  e 


Sgref(kp,oc,ru)  =  (-]  kp  HaiZ,(kpru)  e 
\cj 


(114) 


(115) 


(116) 


(117) 
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Technically,  the  operations  in  (1 14)-(1 15)  represent  the  discussed  deconvolution  in  the 
spectral  domain  with  the  system  functions  defined  in  (1 16)-(1 17).  This  recovery,  however,  is 
not  practical  due  to  the  discussed  (Section  5.1)  finite  support  band  of  reference 

(^rcf^KK/) and  received  (Sft,c,S^ )  signals.  A  practical  reconstruction  can  be  achieved  by 
using  a  matched  filtering  scheme  such  as: 


F/i „ia)(kU!>a,<o)  =  S'nf  (kp  ,a,ru,  a)Snc  (kli:  ,a,ru,  at)  (1 18) 

^Ja\kp,a)  =  SlKj{kp,a,ru)S!l{a,rll,(og)  (119) 


where  (S*er,S*ri,j)  are  the  complex  conjugates  of  reference  functions. 
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6.  Coherent  radar  cross  section  (RCS)  analysis 


An  attractive  feature  of  the  present  approach  is  the  possibility  of  direct  coherent  RCS 
derivations  from  the  knowledge  of  target  function  F(r)  (or  F,  (r)  ).  The  information 
contained  in  F(r)  represents  the  electromagnetic  identity  of  the  target  and,  therefore,  is 
capable  of  producing  accurate  estimate  of  scattering-related  observables  such  as  the  coherent 
RCS.  Furthermore,  derivations  of  the  described  observables  are  not  limited  to  the  far  field  of 
the  scatterer  (or  radiator)  and  may  be  carried  out  in  the  near  field  region  by  using  (1 1)  (or  a 
proper  approximation  of  (1 1))  instead  of  (12)-(13).  Prior  to  such  near  field  considerations, 
however,  a  qualified  description  and/or  definition  of  RCS  in  the  near  field  has  to  be 
established. 


The  bistatic  scattering  cross  section  for  an  incident  plane-wave  E'  can  be  expressed  as  [28]: 


a{6,tj>,9'  ,<j>')  =  lim4^r2 


|E"W)|2 

I  E'(0,0)  | 2 


(120) 


According  to  the  far  field  formulation  (36),  we  have: 


cr  {9,  <p,9‘  ,</>')  =  lim4;ro2 


^  f  f  f  Fy.  (r)  <T/A|u-r|  dr 

Am  JJJ 


kV 

An 


JJJ Fy.  (r)  e-jk]u~A  t/r|2 


(121) 


Applying  the  far  field  approximation: 


f .  u 

\ 

A 

|«w-u-r 

u  =  - 

V  u 

/ 

(122) 


and  coordinate  transformation  from  Cartesian  to  Spherical  systems,  we  find: 
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<t>,  0 '  ,<f> ' )  =  a  e  (0,  <f>)  +  a*  (0,  </>) 


(123) 


where. 


£2  2 

°*  (0’0)  =  J/J  (fT,  cos(^)  cos(</>)  +  F.,.  cos(0 )  sin(^)  -  Fn  sin(0))eJk “r  dxdydz  2 

£2  2 

=  (-F7,sm(</>)  +  FTycos(0))eMr  dxdydz  2 


(124) 


(125) 


In  equations  (124)-(125): 


Fr  (r)  =  Fr  (*>  T*  *)  =  ^  (*,  z)  x  +  FTy  (x,  y,z)y  +  Fl:  (x,  y,  z)  z 


(126) 


and 


e7*u  r  =  exp(yA:  (x  sin(#)  cos(^)  +  ^sin(6>)sin(^)  +  zcos(<9))) 


(127) 


The  integrations  in  (124)-(125)  can  be  analytically  or  numerically  carried  out.  For  a  MoM 
formulation,  the  integrals  are  carried  out  for  an  assumed  set  of  basis  functions  (limited 
number  of  analytical  integrations).  The  RCS  is  then  computed  by  summing  over  the  above 
kernel  moments  using  the  weighting  coefficients  determined  through  MoM. 

It  is  theoretically  possible  to  transform  the  Green’s  function  formulation  (e.g.,  (13))  to  the 
spectral  domain  and  derive  the  RCS  from  the  Fourier  knowledge  of  the  target  function.  This 
would  require  an  equivalent  RCS  definition  in  the  spectral  domain. 
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7.  Conclusion 


Utilization  of  radar  or  SAR  polarimetry  (PolSAR)  for  target  recognition  can  significantly 
enhance  target  detection,  classification  and  recognition  capabilities.  To  make  the  efficient  use 
of  such  techniques  for  remote  sensing,  appropriate  data  exploitation  algorithms  should  be 
employed.  Considering  the  nature  of  a  polarimetric  approach  that  essentially  is  the 
characterization  of  a  target  or  scatterer  through  interaction  with  the  polarized  or  vectorial 
electromagnetic  (EM)  field/wave,  the  required  modeling  analysis  has  to  be  based  on  an 
electromagnetic  theory.  On  the  other  hand,  an  EM-based  approach  to  formulate  the  vector 
field  scattering  is  by  nature  involved  and  complicated.  Noticing  the  need  for  an  optimal 
approach  (i.e.,  rigorous  and  feasible),  the  formulation  here  is  structured  to  satisfy  the  accuracy 
requirements  by  including  the  key  EM  scattering  components  and  developing  a  precise  system 
kernel.  The  feasibility  of  associated  algorithm  implementation  is  addressed  by  formulating  a 
generic  EM  scattering  problem  with  known  (within  other  EM  areas)  analytical  or  numerical 
solutions. 

In  the  present  theoretical  framework,  the  target  scattering  and  classification  problem  is 
approached  by  a  method  which  is  new  within  the  realm  of  conventional  PolSAR  target 
characterization  and  classification  techniques.  The  technique  which  is  based  on  the 
characterization  of  target  scatterers  by  current  vector  functions  (i.e.,  the  principle  of 
equivalent  sources)  and  formulation  of  scattering  using  accurate  Green’s  or  system  kernel 
functions  in  the  spectral  domain,  is  well  founded  in  the  field  of  classical  electromagnetics. 
Introducing  appropriate  basis  sets  representing  various  scattering  features  for  target  function 
decomposition,  characteristics  of  a  target  scatterer  may  be  explored  from  different  aspects  by 
making  use  of  the  derived  accurate  system  kernel.  Theoretically,  there  is  no  restriction  on  the 
nature  and  number  of  features  that  can  be  exploited  as  long  as  they  are  modelled  as  basis 
functions.  Strong  and  insightful  techniques  from  other  domains  of  computational 
electromagnetics,  such  as  microwaves  and  millimetre-wave  structures,  may  be  adopted  and 
tailored  for  appropriate  feature  determination,  analytic  expression  and  effective  target 
function  projection  into  the  feature  space.  Examination  of  the  current  distribution  behaviour 
on  the  target  scatterer’s  surface,  which  is  one  of  the  outputs  of  current  analysis,  can  provide 
good  understanding  of  the  scatterer  physical  characteristics  (e.g.,  shape).  The  current 
distribution  on  the  scatterer  may  be  readily  applied  to  derive  near  and  far  field  characteristics 
and  observable  quantities  such  as  RCS.  The  latter  can  also  be  used  to  extract  target 
distribution  information  from  RCS  knowledge  through  inverse  analysis  and  incorporation  of  a 
priori  information.  Accordingly,  one  may  effectively  utilize  conventional  radar  measurements 
to  validate  and  tune  the  present  scattering  model. 

Although  fundamental,  the  analysis  of  scattering  developed  here  for  polarimetric  radar 
applications  is  general  and  should  be  tailored  for  the  SAR  scenario  of  interest.  The  latter  is 
essentially  done  by  defining  the  appropriate  analytic  basis  functions  that  represent  target 
features  based  on  the  nature  of  scattering  and  the  scatterer.  Decomposition  and/or  projection 
of  the  target  scatterer’s  current  function  into  the  features  (i.e.,  basis  functions)  space  should 
also  be  optimized  by  accounting  for  the  analytic  nature  of  target  and  basis  functions.  The 
future  work  will  be  focused  on  such  target  scattering  feature  classifications  and  required 
analytic  techniques  for  associated  target  function  projections.  A  variety  of  EM  computational 
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techniques  will  be  considered,  tailored  or  further  developed  for  desired  applications.  Results 
of  the  available  trials  and  experiments  will  be  used  and  (if  required)  new  experimental 
schemes  will  be  designed  to  validate  and  tune  the  model  by  taking  advantage  of  the 
established  relations  that  connect  measured  quantities  (e.g.,  RCS)  to  scattering  functions. 
Experimental  data  and  EM  simulations  for  known  target  scatterers  will  be  employed  to 
investigate,  derive  and  designate  analytic  characteristics  to  scattering  features  and/or  classes 
of  features. 
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9.  Annexes 


9.1  Annex  A:  Reference  frame  unit  vector  transformation 


It  is  often  useful  to  express  the  scattered  field  vectors  or  matrices  defined  in  the  radar 
reference  frame  in  terms  of  vectors  in  a  universal  reference  frame,  e.g.,  multi-incident-angle 
analysis. 

A  ^  A 

Assume  (h,  v,  1)  as  the  basis  set  for  radar  reference  vectors.  With  respect  to  a  universal 
reference  frame  (x,y,z) ,  we  have: 


(■ Al ) 


where, 


Uhti  = 

V 

,  U«*  = 

X 

y 

L'J 

(A2) 


f  = 


<cos(0J)cos(0) 

-sin(^) 

vsin(<9()cos(^) 


cos(6>,. )  sin(jzi)  -  sin(#,  )> 
cos(^)  0 

sin(£? )  sin(^)  cos(0,. )  J 


(A3) 


In  (A3),  Oj  denotes  the  incident  angle  and  tf>  represents  the  polar  angle  in  the  transversal 
plane  (to  LOS). 

Similarly,  one  finds  for  the  scattering  matrix: 
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9.2  Annex  B:  Principle  of  stationary  phase 


The  principle  of  stationary  phase  (PSP)  states  that  [1],[34],  for  problems  involving 
integrations  of  complicated  waveforms  with  widely-varying  phase  functions,  the  phase  is 
variational  around  a  certain  stationary  point.  Considering  such  integrand  as: 


Integrand  =  /  (v)  exp(yO(v)) 


(51) 


where  /(v)  is  assumed  to  be  a  slowly-varying  envelope,  the  PSP  is  formulated  as: 


<50(v)  =  0  (52) 


Relation  (B2)  provides  a  powerful  tool  that  can  be  effectively  applied  for  integrations 
associated  with  spectral  analysis  of  complex  functions.  Such  complexity  can  be  gauged  by  the 
time-bandwidth  product  (TBP)  of  the  waveform. 

In  (B2),  the  stationary  point  V  is  considered  to  be  a  vector  in  general.  For  single-dimensional 
problems  (normally  the  case),  the  stationary  point  satisfies: 


ov 


Expanding  0(1/)  in  the  vicinity  of  v0  : 
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(54) 


<D(v)  *  &(v0 )  +  O>0 ) 


(a,>o)s_<p(^)|_ii) 


and  using  the  known  result, 


\eJa{v~v “)'dv  = 


■  \'n 

2IA 

K  «  J 


(55) 


one  finds: 


ff(v)exp(j&(v))ctv 


r  r 


/K)exp(y(D(v/0)) 


(56) 


The  phase  function  is  assumed  to  be  analytic  (near  v0 )  in  the  above  derivations.  The  accuracy 

of  (B6)  depends  on  the  significance  of  the  residual  terms  in  (B4).  Effects  of  the  residuals  are 
inversely  proportional  to  TBP  of  the  waveform. 
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List  of  acronyms  and  initialisms 


CTD 

Coherent  Target  Detection 

DND 

Department  of  National  Defence 

EM 

Electromagnetics 

LOS 

Line  Of  Sight 

MoM 

Method  of  Moments 

PCTD 

Partially  Coherent  Target  Detection 

p.e.c 

perfect  electric  conductor 

p.m.c 

perfect  magnetic  conductor 

PM 

Phase  Modulated 

PSP 

Principle  of  Stationary  Phase 

RCS 

Radar  Cross  Section 

SAR 

Synthetic  Aperture  Radar 

TBP 

Time  Bandwidth  Product 
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